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Abstract 

Yang-Mills theory with flavor quarks in the CIS4 is studied through the dual 
supergravity in the AdSs x S 5 background with non-trivial dilaton and axion. 
The flavor quarks are introduced by embedding a probe D7 brane. We find that 
the dynamical properties of YM theory in the dS4 are similar to the case of the 
finite temperature theory given by the 5d AdS-Schwarzschild background. In 
the case of dS4, however, contrary to the finite temperature case, the gauge field 
condensate plays an important role on the dynamical properties of quarks. We 
also give the quark-antiquark potential and meson spectra to find possible quark- 
bound states. And we arrive at the conclusion that, while the quarks are not 
confined in the dS4, we could find stable meson states at very small cosmological 
constant as expected in the present universe. But there would be no hadrons at 
early universe as in the inflation era. 
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1 Introduction 



Recently, there have been many approaches to QCD from the gauge/gravity correspon- 
dence based on the superstring theory pQ. Especially, the idea, proposed by Karch and 
Katz [2], to add light flavor quarks by embedding D7 brane(s) as probe has stimulated 
many authors who have developed and extended this idea to various lOd gravity solu- 
tions corresponding to the various gauge theories. And many interesting results have 
been obtained for the properties of quarks and and their bound states, mesons, in the 
context of the holography [SlUSISHElEJIin]. 

This approach has been developed in many directions, but the analyses are re- 
stricted to the gauge theory in 4d Minkowski space-time and to the finite temperature 
theory. On the other hand, some holographic approaches to the gauge theory in the 
4d de-Sitter space (dS^ are seen [III [I2l DSl EE] , but it is not enough and we will need 
more study. This situation of the positive cosmological constant can be set for the 
early inflation universe or for the present small acceleration which has been observed 
recently in our universe. From this cosmological viewpoint, it would be important to 
make clear the non-perturbative properties of the gauge theories at finite cosmologi- 
cal constant. It would be a difficult problem to see the non-perturbative behaviors of 
field theories in a curved space from the standard field theory, so this problem is a 
challenging one. 

Here we study this problem from the holographic approach which has been useful 
in the finite temperature case [9]. The bulk solution corresponding to the gauge theory 
in dS4 is obtained from type IIB string theory with dilaton and axion under the five 
form flux. And the D7 brane is embedded in this background as a probe to introduce 
the flavor quarks. In the bulk background, there appears a horizon as in the finite 
temperature case. And we find an attractive force between D3 and D7 branes, then the 
chiral symmetry is preserved as in the high temperature phase. But a phase transition 
observed in the high temperature case is seen only for the case of large gauge field 
condensate. For small gauge condensate, it is deformed to a different situation which 
has not seen in the finite temperature case. 

In order to study the quark confinement, we obtained the potential between quark 
and anti-quark through the estimation of the Wilson-Polyakov loop. We find a max- 
imum distance between quark and antiquark to maintain the U-shaped string state. 
Above this length, the quark and antiquark are separated as independent particles. 
And the energy of each single quark (or anti-quark) state is finite. This point is as- 
sured through the estimation of the effective quark mass. Then, in this sense, we can 
say that the theory is in the deconfinement phase. We compare these results with the 
one given at the finite temperature. 

Furthermore, we estimate the meson masses through the fluctuation of the D7 
brane, and we find interesting spectra for light mesons. However, all the states would 
disappear at large cosmological constant since they becomes unstable and decay to free 
quarks and ant-quarks. A similar phenomenon for the baryon spectra is also seen by 
studying the energy of the D5 baryon wrapped on S 5 which is regarded as baryon. 

In section 2, we give the setting of our model for our study, and a problem of the 
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phase transition is discussed by solving the embedding of the D7 brane. In section 
3, effective mass of the quark and the quark- antiquark potential are studied through 
the Wilson Polyakov loop estimations. In section 4, the possible bound state for the 
mesons and baryons are discussed. The summary is given in the final section. 



2 Background geometry and D7 brane embedding 

We solve the equations of motion for lOd IIB model retaining the dilaton $, axion x and 
selfdual five form field strength F( 5 ). Under the Freund- Rubin ansatz for F^, -F Ml ... M5 = 
-\fK/2 e Atl ... M5 [151 US], and for the lOd metric as M 5 x S 5 or ds 2 = g M N<lx M dx N + 
gijdx l dx\ we can find the solution. The five dimensional M5 part of the solution is 
obtained by solving the following 5d reduced action, 

S = ^ / d 5 x^ (R + 3A - l -{d§Y + \e 2 *{d X ) 2 

which is written in the string frame and taking a' — g s — 1. Then, we obtain, (see the 
Appendix A) 

ds 2 Q = GcMjydX^ dX^ 




A 2 [-dt 2 + a(t) 2 (dx 1 ) 2 ) + ^dr 2 + R 2 dnl\ , (2) 



q 1 - (r /r) 2 /3 



r 



X = -e + Xo , (3) 



r 4 (1 - (r /r) 2 ) 3 

A = l-(^) a , a {t)=e 2 % t M) 
r 

where M, iV = ~ 9 and R = v^A/2 = (AnN) 1 ^. As for the integration constants, r 
denotes the horizon point, and it is related to the 4d cosmological constant A as 

A = 4£. (5, 

And q is a constant which corresponds to the VEV of gauge fields condensate |9]. And 
other field configurations are set to be zero here. We notice that there is no singularity 
at the horizon r in this case. This situation is similar to the black hole configuration 
used for the finite temperature gauge theory. 

In the present configuration, the four dimensional boundary represents the infla- 
tional universe characterized by the 4d cosmological constant A. So the above bulk- 
solution describes the gauge theory in the inflation universe. Firstly, we examine the 
properties of the flavor quarks in this gauge theory, and it is obtained by solving the 
embedding problem of the D7 brane in this background. 

D7 brane embedding and phase transition: The D7 brane is embedded as 
follows. The extra six dimensional part of the above metric ([2]) is rewritten as, 

^dr 2 + R 2 dn 2 5 = ^ (dp 2 + p 2 dn 2 3 + (dX 8 ) 2 + (dX 9 ) 2 ) , (6) 
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where r 2 



(X 8 ) 2 + (X 9 ) 2 . And we obtain the induced metric for D7 brane, 



e" ,2 \^- 2 A 2 (-dt 2 + a{t) 2 {dx*) 2 ) + 



R 



— ((i + (d p w) 2 )dp 2 + P 2 dn 



(7) 



where we set as X 9 = and X 8 = w(p) without loss of generality due to the rotational 
invariance in X 8 — X 9 plane. Then the embedding problem is reduced to obtain the 
solution for the profile function w(p), and it is performed as follows. 

The brane action for the D7-probe is given as 

Sl)7 



-T7 / cm e 



det (G ab + 2na'F ah ) - ^e 11 "" 8 ^. 



-is 



+ M^p [C (o ]APAF 



where F ab = d a A b - d b A a . Q ah = d^X M d^X N G M N (a, b = ~ 7) and r 7 = 
[(2ir) 7 g s a' 4 ] _1 represent the induced metric and the tension of D7 brane respectively. 
And PfC^ 4 -*] denotes the pullback of a bulk four form potential, 



C (4) 



i? 4 



dx° A dx 1 A dx 2 A dx 3 



(9) 



The eight form potential A^...^, which is the Hodge dual to the axion, couples to the 
D7 brane minimally. In terms of the Hodge dual field strength, Frg\ = dA^ [17], the 
potential A( 8 \ is obtained. 

When the gauge potentials A a are neglected, the action is abbreviated as 



S D7 = -r 7 J d^ (e- 
and by taking the canonical gauge, we arrive at the following D7 brane action, 
Sbt = ~r 7 J d 8 ^p 3 a(t) 3 (A*e*y/l + (W) 2 - C 8 ) , 

C 8 (r) = f dr' A\r')d r , (exp($(r / ))) = \. 
Then, finally, the equation of motion for w(p) is obtained as, 

w 



(10) 



(12) 



p + w w' 



$' - a/1 + (w') 2 ($ + 41ogA) / 



1 



1 + (w 



r\2 



w' f - + ($ + 4 log A)' J + 



l + (w 



A2 



(13) 



where prime denotes the derivative with respect to p. By solving this equation we 
find the profile of the embedded D7 brane and then we find simultaneously the quark 
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Fig. 1: Potential of to, V(w) = A 4 e* + C 8 , for g = 0, 0.3 and 10 with R 
r = 1.0. 



properties, the quark mass m q and the chiral condensate where $ denotes the 

quark field. This point is shown below through various explicit solutions. 

Trivial Solution and Chiral Symmetry: Firstly we consider the asymptotic 
solution of w for large p. It is obtained usually as the following form 



w(p) ~ m q + — , 
P 



(14) 



where m q and c are interpreted from the gauge/gravity correspondence as the current 
quark mass and the chiral condensate, respectively, since the field <ft 8 is corresponding 
to the conformal dimension three operators of the gauge theory. General solution of 
([TBI is characterized by these two arbitrary parameters, m q and c, since Eq. f lTBl is the 
second order differential equation of w although it is highly non-linear. However, for a 
fixed m q , we find one solution which is meaningful from the point of view of holography 
[9]. In other words, c is determined when m q is fixed. This can be interpreted in the 
gauge theory as that — c = is determined by a theory with a parameter m q , the 

quark mass. Then the solutions for w are characterized only by the quark mass m q , 
and the vev of quark condensate is determined uniquely by m q . 



However we must be careful in using the above asymptotic form (jl4"|) in solving the 
equation (fl3j) since the form (fT4l) is useful only for the case that CFT is realized on 
the boundary r ~ p — ► oo. In the present case, the geometry of the 4d boundary is 
dS4 and the conformal symmetry is broken there, then the form (|14p is not useful in 
getting the solutions of f|T3|) . 

Actually, we find a solution of m q = 0, and arbitrary c when we solve the Eq. fflBl 
by expanding it in terms of the power series of 1/p 2 with the asymptotic form of (JHJ). 
Then, we find that the meaningful solution of m q = is nothing but the trivial solution, 
w = 0. The result is however non-trivial and important, because it implies that the 
chiral symmetry is preserved in the dS4. 
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This point is understood from the fact that the attractive force is working between 
the D3-branes at the horizon and the D7 brane at X 8 = w. The force between them is 
obtained from the potential of w, which is obtained from the D7 action (iTTj) by setting 
w' = and remembering r 2 = p 2 + w 2 as follows, 

V(w) = r 7 (A 4 - C s ) . (15) 

It represents the potential of w at fixed p, and typical examples for different qs are 
shown in the Fig. [TJ 

Then the force F between D3 and D7 branes is obtained 

F=-%- = -Sr 7 ^°A*e* < 0, (16) 

and we can see that the force is attractive at any point of p. 

Then we can understand that the c must be negative for any solution of w. As a 
result, in the case of m q = 0, we arrive at the trivial solution, m q = c = or w = for 
non-zero A. This situation is similar to the case of finite temperature gauge theory, but 
we find several different features through the non-trivial solutions between the finite 
temperature and finite cosmological constant case. 

Non-trivial Solutions and Phase Transition: However, the trivial solution is 
a part of solutions. In fact, we could find many non-trivial numerical solutions, whose 
behavior seems to be consistent with the asymptotic form of (|14p at a glance. Some 
of them are shown in the Fig. [2]|3] for different qs. In order to obtain these solutions 
in terms of the power series expansions, we must add correction terms of log(p) to 
the chiral condensate c and to the other coefficients of the power series of l/p 2 . This 
improvement is naturally expected since the VEV, c = —(^^, would receive this 
type of loop corrections due to the conformal non-invariance of the gauge theory in the 
present background dS4. 

Then, instead of (j!4p . we find the following asymptotic form 

/ \ , c - Am q rl log(p) 

w(p) ~ m q + \ , 17 

p 2 

where m q and cq are arbitrary at this stage. However, cq would be determined by m q 
and other parameters of the theory for the meaningful solution as stated above. And 
the solutions of w(p) and c for various m q are shown by separating the solutions to 
three groups depending on the value of q. 

As for the c, in the present case, it depends on log(p) and diverges at p = oo. So we 
need an appropriate subtraction or a renormalization for this quantity [TJ]. However, 
we estimate it as 

c = c - 4m,r 2 log(p cutofr ) (18) 
here by introducing an appropriate cutoff p C utoff f° r the simplicity. 

In any case, there is a horizon at p 2 + w(p) 2 = r 2 in the p — w plane, and the 
embedded solutions are separated to two categories by the infrared end point of w(p). 
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Fig. 2: Embedding solutions w and the chiral condensate c for q = and tq = R = 
1.0. For c, the cutoff is taken at p = 100, and the dashed line shows the transition 
point from group (a) to (b) at m q = 2.639. 




Fig. 3: For q = 0.3 and other parameter settings are the same with the above Fig.[2j 
The dashed lines in the right figure show the two critical points, m q = 2.683 and 

= 771, = 2.118. 
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Fig. 4: For g = 10 and other parameter settings are the same with the above Fig.[2j 



In the first group, it is above the horizon (group (a)), and for the second, it is on the 
horizon (group (b)). This situation is also seen in the case of the finite temperature 
gauge theory. 

The solutions are shown in the Fig. [20 for different qs and we find that all the 
solutions of w decreases with decreasing p, which means c (= Co — 4m g rQlog(p)) < 0, 
for any solution. This is, as mentioned above, due to the attractive force between D3 
and D7 branes, and c = for m q = as expected. In the Fig. CEJ we show the potential 
for the three different q cases for our understanding of this point. In other words, 
the chiral symmetry is preserved for m q = in any case. The numerical results are 
summarized as follows. 

(i) q = 0: Firstly we consider the case of q = or zero gauge field condensate, 

shown in the Fig. [2, (Fjlv) = 0- m this s t a te> the quarks are not confined even in the 
limit of A = 0. The second point to be noticed is that we could not find any jump of 
solutions when it changes from the group (a) to the group (b), and c changes smoothly 
at the transition point of m q . 

For finite q, the solutions for the group (b) are strongly restricted. Namely, the end 
point value of wo = w(po), where po denote the value of p on the horizon, is constrained 
as (see Appendix B) 

Wo < 7!o (19) 



and the maximum of Wq is realized in the limit q — > but q ^ 0. We find solutions 
of m q < 0.96 (in Gev) for q = 10~ 7 and r = 1. On the other hand, the solution for 
q = and r = 1 has m q > 2.4, so there is a gap for the quark mass of groups (a) and 
(b). This gap decreases with increasing q, and we summarize the numerical results for 
q > as follows. 

(ii) Small g(< 2): We show in the Fig. [3] the solutions for non-zero but small q, 
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namely at q = 0.3. The p dependence of the solutions are similar to the q = case, 
but there is a forbidden region of solutions for 2.118 < m q < 2.683 (in GeV) between 
the two solution groups of (a) and (b). This phenomena is characteristic to the present 
case of dS 4 , and it has not been seen in the finite temperature case. This gap gradually 
decreses with increasing q and disappears for q > 2 as shown below. 

(iii) Large q(> 2): In the Fig. HJ we show the solutions for large q{= 10) case. In 
this case, the behavior of the solutions is similar to the case of the finite temperature 
gauge theory. Namely, the forbidden band of the quark mass or a gap of tu(oo) at 
the boundary seen in (ii) case disappears, but the gap at the infrared end point of w 
still remains and a phase transition as seen in the finite temperature case is expected 
when the solution translates from group (a) to the group (b). This is seen as a gap 
of the infrared end point w for the solution of common m q at the transition point. In 
the present example, it is realized at m q = 2.94(= m c ). This transition is also seen 
through the chiral condensate c, which is actually seen from the Fig. H] as a jump of c 
at m q = m c . In the figure showing c, the upper one shows the one of solution group 
(a) and it begins with a gap at m q ~ 2.9. And there is a narrow double valued region 
near the point A. In the multiple valued region of c, we choose the value of c, whose 
D7 energy is the minimum, as the plausible solution. The lower curve is obtained from 
the solutions of group (b) and it ends at m q ~ 3.2, the point B. The similar gap is seen 
also in the mass spectra of mesons as seen below. 

Next, we study the reguralized D7 energy E D7 defined as 

S D7 = -T 7 J d 7 ^a(t) 3 E m (20) 
E D7 = r dp p 3 fv4Vyi + (w 1 ) 2 - Cs) , (21) 

in order to see the order of the phase transition which is clearly expected in the above 
calculation of case (iii) for g(> 2). The regularization for the D7 energy for the finite m q 
solution is usually performed by subtructing the one of m q = [3] , and the regularized 
energy E^f(m q ) for the mass m q is given as 

E%*(m q ) = E m (m q ) - E D7 (0) (22) 




where p mu tch represents the position where w is approximated well by its asymptotic 
form ( FlTl) . In the present case, however, we need an extra subtraction term due to 
the cosmological constant tq. In fact, in the large p region, the integrand ALd7 is 
expanded by using asymptotic form of w(p) as follows 

2 2 

AL D7 = 4^ + 0(p- 3 ) (24) 
P 
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Fig. 5: The regularized energy for q = 0, and other parameter settings are the 
same with the above Fig. [2j 



e r 




Fig. 6: The regularized energy for q = 0.3, and other parameter settings are the 
same with the above Fig. [2j 



Er 




Fig. 7: The regularized energy for q — 10, and other parameter settings are the 
same with the above Fig. [21 The points A and B are the same points with the 
onein Fig. |H 
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so we subtract here the first term in the estimation of E^!f(m q ). Then we obtain a 
small correction from the third term of Eq. (1231) . 

The numerical results are shown in the Fig [SHU In the case of q = 0, we can see 
a small cusp-like structure at the transition point between the two types of solutions 
(a) and (b). This might imply some kind of phase transition. On the other hand, for 
q = 10, we find definite energy gap between the solutions of two groups at the same 
value of m q in their overlap region. We then expect a first order phase transition from 
a solution of group (b) to the one of group (a) at an appropriate point of m q . The 
details of these transitions observed here will be given in the future. 



3 Quark- ant iquark potential and quark confinement 

We study a gravity description of quark- ant iquark potentials in order to study the quark 
confinement. Before giving the calculation, we briefly review how quark-antiquark 
potentials described in the context of the gauge/gravity correspondence. 

Consider the Wilson-Polyakov loop, W = (1/JV)'IWV/ A ° d \ in SU(N) gauge the- 
ory, then the quark-antiquark potential V qq is derived from the expectation value of a 

parallel Wilson-Polyakov loop as (W) ~ e~ Vqq ^ ' dt . Meanwhile, from the gravity side 
the expectation value is represented as 

(W) ~ e" 5 , (25) 

in terms of the Nambu-Goto action 

S = J drda^J -det h ab , (26) 

with the induced metric h a b = G lxv d a X ll di ) X 1 ' ', where X^(t, a) denotes the string coor- 
dinate. 

Then the quark-antiquark potential can be calculated by setting the configurations 
of string coordinates under the background geometry given above. To fix the static 
string configurations of a pair of quark and anti-quark, we choose X° —t — r and de- 
compose the other nine string coordinates into components parallel and perpendicular 
to the D3-branes: 

X = (X,|, r,rfi 5 ). (27) 
Then the Nambu-Goto Lagrangian in the background ([3]) becomes 



27ra 



— [da e* /2 jA(r)V 2 + r 2 A(r) 2 fl 5 ' 2 + ( -) A(r) 4 a(t) 2 X|,' 2 , 
to! J V \RJ 



(28) 



where the prime denotes a derivative with respect to a. The test string has two possible 
configurations: (i) a pair of parallel string, which connects horizon and the D7 brane, 
and (ii) a U-shaped string whose two end-points are on the D7 brane. 

Effective quark mass and confinement: Firstly consider the configuration (i) 
of parallel two strings, which have no correlation each other. The total energy is then 
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two times of one effective quark mass, m q . As mentioned above, it is given by a string 
configuration which stretches between r and the maximum r max , so we can take as 

r = cr, X|| = constant, Q§ = constant. (29) 

Then rh q is obtained by substituting (|29j) into (1281) as follows, 

E = / dr e* /2 A(r) = 2m q , (30) 

Tia' Jr 

where r rmax denotes the position of the D7 brane. 




Fig. 8: 2m q are shown for R — 1 (GeV 1 ),r rnax = 3 and a' = 1 (GeV 2 ). The three 
curves are corresponding to the case of q=0, 0.3 and 3 respectively. 



The integrand e $//2 A(r) is diverges at r = r as 1/y/r — r , but we find that the 
contribution of this part to the integration vanishes, 

/ ^7^== = 2v / ^a o =0. (31) 

Then we find rh q < oo for finite r or A. This means that the quark is not confined 
in this case since single quark could exist. On the other hand, we find fh q diverges for 
r = then the quark confinement. 

While, in the above discussion, q ^ is essential, we find for q = 

2m q = _l_ ( r ^x-r ) 2 (32) 

TTOi '"max 

Then the quark is not confined in this case even if A = 0. 
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As for the cosmo logical constant dependence on m q , the numerical calculations are 
shown in Fig. [HJ From this Fig., we find that (i) fh q increases with q at any point of A, 
and (ii) it decreases with increasing A for any q. 

U-shaped Wilson-Loop: We now turn to the U-shaped configuration, 

X|| = (<t,0,0), Q 5 = constant. (33) 

The equation of motion derived from the Lagrangian (1281) with the configuration fl33l) 
are solved by 

e* /2 , 1 (^) A 3 (r)a 2 (t) = constant. (34) 

^f{r/R) 4 A 2 (r)a 2 {t) + {dr/daf \R' 

The midpoint r min of the string is determined by dr / da\ r=Tmin = 0. Then the distance 
and the total energy of the quark and anti-quark are given by 

L = 2R / dr (35) 

Jr min r 2 A(r)a(t) v /e*Wr 4 A(r) 4 / {e^ r ^h^ in A(r min ) 4 ) - 1 

E = -L r~ A{r)e ° M ' 2 (36) 

V/1 - e*< r "->r 4 min A(r min )*/ (^Wfr)') 



TTOt' Jr. 



Here we study the time independent distance L defined as L = aL instead of L given 
above. The numerical results are shown in the Fig. [9] 




Fig. 9: Plots of E vs L at g = (the left figure) and = 0.3 (the right figure) (GeV 4 ) 
for R = l(GeV" 1 ), r max = 3(GeV" 1 ) and a' = l(GeV" 2 ). The solid and dashed curves 
represent the case of A = and A = 4 (GeV), respectively. The vertical solid and dashed lines 
represent the energy of two parallel straight strings. 



Two figures show the dependence of the energy E on the distance L at the selected 
cosmological constant A and q. For A = 0, the well-known results are seen for q = 
and finite q; (i) For q = 0, £7 oc 1/L at large L and E oc m 2 L at small L [3]. And for 
finite q [9], E oc ^fqL at large L and E oc m 2 L at small L. The behaviors at small L 
are the common since the same AdS limit is realized there. 
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For finite A, there is a maximum value of L (= L max ) to find the U-shaped con- 
figuration. Namely, the U-shaped configuration disappears for L > L (= L max ). The 
similar behavior is seen also in the case of the finite temperature [21]. In this sense, 
the theory in dS4 is in the quark deconfinement phase as in the finite temperature case. 
However, we notice the following difference. In the case of finite temperature, there are 
two possible U-shaped string configurations at the same values of L (< L max ), but in 
the case of finite cosmo logical constant, U-shaped string configuration is unique at a 
given value of L(L < L max ). And at L — L max , the energy of this string configuration 
arrives at 2fh q . Then, this implies that the U-shaped string configuration is broken for 
L > L max to decay to free quark and anti-quark. On the other hand, an unstable U- 
shaped string configuration is allowed for the finite temperature case even if E > 2rh q , 
and, just in this energy region, the other U-shaped string configuration is formed. 



In our model, L max is given as, 



L max = lim L , (37) 



and for q = 0, and we obtain 



Hm / dx- jL-l) 3^ 

y/\Xo-*lJx h (1 - X 2 ) \X 4 (1 - X^Y 

7T 1 



2^ 



(39) 



where 



X — , X a — , X\j — . J 



While the value of L max in general depends on the lower bound Xb, however the ap- 
proximate value seems to be independent of it. This is because of the fact that the 
integration is approximated by the value near the upper bound, x a , of the integration. 
And we can show the x^ independence of L max directly by 

OXb 

The numerical result is consistent with the one given in 
As for q ^ 0, we get 



-X 2 ) 3 - 

^\ C-l J Xb ^ (1 - X 2 ) V X 4 (1 + ^4|zip) (l _ X 2)4 

Vx ' 

where q = qjr\. In this case also, we can show the r max independence of L max by (HTj) . 
but the q independence of the above result is non-trivial. This point is understood from 
the fact that the dominant part of the integrant in this integration is independent of q. 



(43) 
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As a result, we get the result which is independent of q and r max . But the approximate 
value is about four times of the one with q = 0. This is reduced to the difference of 
the quark- antiquark potentials at A = for q = and q ^ 0. For the latter case q ^ 0, 
we obtain the linear potential at A = 0. Their numerical values are also plotted in the 
Fig.dU]forr max = 10. 



(1/7T) V^Lmax 

31 



2.5 
2 - 

1.5 
1 ■ 

0.5- 



q=0 . 3 



q=0 



0.5 1 1.5 2 2.5 3 3.^ 

Fig. 10: Functional relation between ^^ max and A for q = and q ^ 



4 Possible hadron spectrum 

As shown in the previous section, while the theory is in the quark deconfined phase, 
the U-shaped string configuration exists. However, the energy of the configuration is 
bounded from the above by 2m q which is the energy of the parallel string configuration. 
And its value depends on the cosmological constant, A, and q. This situation implies 
the existence of stable meson states, the bound state of quark and anti-quark, when 
their mass is below 2m q . We study this point by calculating the meson mass and 
comparing it with 2fh q . And we find a stable meson when its mass is smaller than 2m q 
since this meson state can be expressed by the U-shaped string configuration, namely 
SIS db bound state of quark and anti-quark. 



4.1 Meson spectra 



The meson spectrum is obtained by solving the equations of motion of the fields on the 
D7 brane. According to [18], firstly we consider the fluctuations of the scalar mesons 
which are defined as, 

x 9 = 9 , x 8 = w( P ) + 4> s . 

And writing the wave functions in the following factorized form, 

4> k = <j ) k (t,x l )<j ) k (p)y l (S 3 ) (k = 9, 8) 
we get the linearlized field equations for 9 (p) and 8 (p) as follows 



1 



dy + -^d p {L )d P 4> 9 



[1 + w 



I 2\ 



R ^ 77lg 



1(1 + 2) 
P 2 



2K (l) 
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+ (i WY^/ 



r dr 



and 



^1 



" i^ 4 m 2 Z(Z + 2) 



eM 4 



d^A 4 ) 



(44) 
(45) 



r 7 A 2 



2(1 + w' 2 )(if (1) + 2w 2 AT (2) ) 



+(l + w 



I 2\3/2 



9$ <9 2 $\ w 2 d$p 2 



Li 



= -2—d p (L w w'K (1) )<f) 8 



(46) 
(47) 



Where four dimensional mass mg and m§ are defined by 

^ k (t,x i ) + 3-^ k (t,x i ) + ^<l> k (t,^) = -mU k (t,x i ). (k — 9, 8) (48) 
a ar 

In deriving the above equations of <f) 8 and 9 , we used 

r 2 = p 2 + (0 8 ) 2 + ((f) 9 ) 2 + w 2 + 2w<f . (49) 

But we should notice here that the variable r in the above field equations is understood 
as r 2 = p 2 + w 2 since we are considering the linearized equations. 



In Figs. 11(a) 11(c) 11(e) , the numerical results of the mass eigenvalues, wig (small 



point/thin line) and m 8 (large point/thick line), are plotted as functions of m q . These 
values are all for the nodeless solutions, i.e. for the lowest mass state. The eigenvalues 
for q 



0, Fig. 11 (a 



show smooth variation with respect to m q as expected from the 
solutions of its profile function (see the Fig. [2D- For a non-zero but small value of q, 
q = 0.3, the mass eigenvalues are shown in the Fig. 11(c) for the region except for the 



forbidden range, 2.11 < m q < 2.68, which has been found and shown in Fig. [3] through 
the solutions of w. In Fig. 11(e), the results for q = 10 are shown. In this case, we 
could find the mass gap at the transition point of m q (shown by dashed line), which is 
found through the solutions of w (see Fig. @J. 

The next problem is to compare the mass eigenvalues wig and m§ obtained with 2rh q 
[ above in order to see the stability of the meson states. In Figs. 11(b) 
m g , m 8 and E = 2m q are plotted with respect to r = y/X/2. In the cas 



of q = (Fig. 11(b)), mg and ms are larger than 2fh q at any point of A. Then there 
is no stable meson state for q = 0. This is reasonable since the quark is not confined 
at A = for q = and positive A works to destroy the confinement phase. On the 
other hand, for q > 0, the quark is confined at A = 0, then we expect stable mesons 
for small A. In fact, in the case of q > 0, we could find the small A region where the 
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I . m 5-8 




(e) (f) 

Fig. 11: mg i8 vs m q for r = 1.0 and R = 1 and (a) q — 0, (c) q — 0.3, (e) g = 10. 
Small points denote mg while large points denote mg ((a)). Thin lines denote mg while 
thick lines denote m$ ((c), (e)). Each vertical dashed line of the right of (a), (c) and 

(e) separates the region which touch the horizon and the region which do not touch the 
horizon. And m 9j8 vs r = a/A/2 for R = 1 and w(0) = 3.0 and (b) q — 0, (d) q — 0.3, 

(f) q = 10. Thin line denotes m 9 while mid-thick lines denotes m 8 . Thick line denote 
E = 2m q . 
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meson masses are smaller than 2m q . From Figs. 11(d) and 11(f) , the stable regions are 
found as r = \f\j2 < 0.1 (0.2) for q = 0.3 (q = 10). The critical value (or the upper 
bound) of A increases with q as expected. 

In our present universe, the value of g 1//4 is expected as the order of hadron mass 
since it provides the QCD string tension [9]. Meanwhile, A is expected to be very small 
from the recent observation of small acceleration even if it exists. Then the mesons 
in the present world would not be affected by this small A and all the mesons are 
stable. On the other hand, there would be no hadrons in the early universe which is 
inflationally expanding under a large A. And the quark and gluons would move almost 
freely in spite of their strong interactions which could combine them to make a bound 
state at A = 0. 

As for the vector field, its equation of motion is obtained from the action (jSJ) as [18J 

d a (^-detg cd F ab ) - Ap ( p2 + w2 h b i%A k = , (50) 

where e 1 ^ is a tensor density (i.e., it takes values ±1). The second term comes from 
the Wess-Zumino part of the action, proportional to the pullback of the RR five-form 
field strength, and is present only if b is one of the S 3 indices. And g a b are metric in 
the Einstein frame. 

Because of the existence of the cosmological constant, Lorentz invariance is broken. 
So, the resulting equations are different for a time component, space components and 
the components with S 3 indices even if A p = gauge is taken. But the final results 
are expected to be not so different from the scalar meson spectra. So, we abbreviate 
the analysis of vectors here. 



4.2 Baryon 

It has been shown that baryons correspond to D5-branes wrapped around the compact 
manifold M5 [TpJ, [20]. Here we assume it to be S 5 . The brane action of such a D5 
probe is 

S D5 = -T 6 / d^Vg , (51) 

where (&) = (X°, A 5 ~ A 9 ), r 5 represents the tension of D5 brane, and Q = — det(^j) 
for the induced metric Qij = d^X M d^X N GMN- The mass of the wrapped D5-brane 
is then 

M D5 (r) = t^VG = r 5 7r 3 i? 4 rA(r)e*/ 2 . (52) 
Before seeing the A dependence of this quantity, we consider the case of A = 0, 

M D5 (r) = r^R 4 r^jl + ^- A . (53) 

This has a global minimum MD5(r m ; n ) = T57r 3 i? 4 (4g) 1//4 , which is regarded as the baryon 
mass, at r = r min = g 1 / 4 .Thus, the baryon mass is also induced by the q, i.e. by the 
gauge-field condensate. This is consistent with the fact that the QCD string tension is 
given by q in the present model. 
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M D5 




1.5 2 2.5 3 

Fig. 12: The D5-brane mass Mn 5 (r) as a function of r. Here we set q = 1, R = 1 
and r 5 = l/7r 3 _R 4 . The curves show the of M D5 | A=4 , M D5 | A=4 . 41 + 1, M D5 | A=4 . 67 + 2, 
Md5|a=5.57 + 3 and Mds|a=7.29 + 4 respectively, from bottom to up. 



Fig. [121 shows the r dependence of Mj) 5 (r) for five values of A. For small side of A, 
there exists a minimum at an appropriate point of r, and it disappears for large enough 
value of A. Thus, there is a transition point of A of confinement for the baryon as seen 
in the meson case. 



5 Summary 

The Yang-Mills theory with light flavor quarks is investigated in the inflationally ex- 
panding 4d space-time or in the time-dependent dS 4 space-time. The flavor quarks 
are introduced by embedding the D7 brane as a probe in the background of the dual 
supergravity. The lOd background is deformed from AdSs x S 5 by the dilaton and 
axion, and its 4d boundary of AdSs is set as the time-dependent dS 4 space-time with 
a 4d cosmological constant which is given as an arbitrary constant parameter in our 
model. 

In this model, the conformal invariance is broken even at the ultraviolet boundary, 
then, in obtaining the D7 embedding, the asymptotic form of the profile function w(p) 
must include logarithmic terms coming from the loop-corrections. Here we find the 
following form 

c - Am q rl log(p) 



w(p) ~ m q + 



P 2 



at the lowest order of large p limit, and this implies the vev of the bilinear of quark 
fields, \I/^> , receives the loop correction proportional to the cosmological constant A 
and the quark mass m q as 

C_ Cp_ 

R 4 ~ R A 



- = _ = _JL _ mq x log(p). (54) 
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This kind of correction would be expected in other quantities also, and we like to turn 
on this point in the future. 

In terms of the above asymptotic form, we obtain various solutions of w(p) with 
different values of m q and Cq. And we find c < for any solution of m q > and c = 
for m q = 0. This implies that the chiral symmetry is kept being unbroken for dS4. The 
solutions are separated to two groups by their infrared end point whether it is above 
the horizon or just on the horizon. And when the solution is switched from the one 
group to the other, three kinds of behaviors are observed. For q — 0, the transition of 
the solutions is smooth, but there is a forbidden region of m q for the case of small q. 
For large q > 0, on the other hand, we find a phase transition, at an appropriate value 
of m q , and this is similar to the one observed in the finite temperature gauge theory. 

In order to see the quark confinement, the Wilson- Polyakov loops are studied. Our 
model for q > shows the quark confinement at A = since we find a linear rising 
potential with respect to the distance between quark and anti-quark, L, in this case. 
On the other hand, in dS4 or for A > 0, the potential increases with L, but it disappears 
for large L, L > L max . And we find that the energy of quark and anti-quark system at 
L = L max is equal to the one of two parallel strings, which connect horizon and the D7 
brane. This means that the quark and anti-quark do not make the bound state any 
more for L > L max and they move freely and independently. In this sense, we can say 
that the gauge theory in dS 4 is in the quark deconfinement phase. 

We should notice that, in the present finite A case, the U-shaped string configuration 
is unique at a given value of L. This point is in sharp contrast to the case of finite 
temperature case, where there are two possible U-shaped string configurations at the 
same values of L. However, the energy of one of them is always higher than the 
one of the two parallel string configuration. In this sense, the stable U-shaped string 
configuration is unique in both cases, i.e. in the finite temperature and finite A cases. 

While the gauge theory of dS 4 is in the deconfinement phase, we expect that some 
meson states are stable for small A due to the following reasons. The value of L max 
proportional to A -1 / 2 , so we will find a U-shaped string with long L at small A. Sec- 
ondly, the energy of the parallel string configuration 2m q increases with decreasing A. 
Actually, rh q becomes infinite at A = for finite q. Then at small A, 2m q exceeds 
the meson mass considered. In order to assure this point, the spectra of mesons are 
examined through the fluctuations of D7 brane. Then, we can show that any meson 
state for a definite quark mass becomes stable when we take A to be small enough. 

As for the baryon, which is identified as D5 brane wrapped on S 5 , its energy is 
obtained by the D5 brane action. We could find that their mass is induced by the 
gauge condensate q in our model. For finite A, in the small side of A, there exists a 
minimum of the D5 energy at an appropriate point of r, and it disappears at large 
value of A. In this sense, we expect a transition point of A of quark confinement for 
the baryon. The details of the transition point will be given in the future work. 
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Appendix A 

We start form the type IIB supergravity with the following bosonic action, 

where other fields are neglected since we need not them. By taking the ansatz for 
F(S), F Ml ... Mfi = e M1 ... M5 jEl HE], and for the lOd metric as M 5 x S 5 or ds 2 = 

gMNdx M dx N + gijdx l dx ] , the equations of motion are given as 

Rmn = ~^9mn + ^d M $d N $ - -b^OmX^nX , (56) 

-t=d M (V=99 MN d N ®) = -e 2 *d MX d NX g MN , (57) 

-t=d M (V^g MN e 2 "d NX ) = , (58) 

V^g v 1 

Rij = Agij ■ (59) 

Using the ansatz, x — — e~*+ constant, the equation ( |56i) is written as 

Rmn = — A^mat- 
Then the metric of M5 part is solved in the following form [22] 

ds 2 {5) = ^A 2 (-dt 2 + a(t) 2 (dx 1 ) 2 ) + ^dr 2 (60) 

A =l-(^) 2 , a (t) = e 2 ^. (61) 
r 

Using this result, $ in the Eqs. fl57l) and fl58l) are solved as given in the section 2. 
We should notice that the integration constants are set appropriately in solving the 
equations. 

Our model is based on type IIB supergravity, and we solved the equations of mo- 
tion with dilaton and axion. In the present case, the model is reduced to 5d gauged 
supergravity. However, the solution given above breaks the supersymmetry since it is 
not a solution of the first order equations, which are written in terms of the superpo- 
tential W / ($, x), which is a constant in the present case, of the 5d theory [23]. This is 
consistent with the fact that there is no supersymmetric theory in dS4. We notice that 
the supersymmetry is also broken by our D7 brane embedding since the k symmetry 
of the D7 brane action is lost [9]. 
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Appendix B 



Here we show the constraint given in Eq.( fT9l) in the Sec. 2. Suppose some solution u>(p) 
of the group (b), then its Infrared end point sits on the horizon at p = po- Then we 
expand w(p) near the end point as 

w(p + e) = w(po) + w'(p )e-\ • (62) 

By substituting this asymptotic form into the equation of w, Eq. (113p . we obtain 



p + 3^/r 2 -9 Wo 2 

w (Po) = r 2 in .„2 > ( 63 ) 

'o ~~ ru«; 

where w = u>(po) and we used p\ = r 2 , — Wq. Since the physical solutions should be 
monotonically increasing with p in the present case, then we obtain wq < r / V4~0. This 
constraint does not depend on the value of q, while this does not appear for q = 0. In 
other words, this behaviour is given by the non-trivial dilaton or the gauge condensate 
in the corresponding 4d gauge theory. 
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